Boundary conditions for coupled quasilinear wave equations with application to isolated systems

We consider the initial-boundary value problem for systems of quasilinear wave equations on domains of the form $[0,T] \times \Sigma$, where $\Sigma$ is a compact manifold with smooth boundaries $\partial\Sigma$. By using an appropriate reduction to a first order symmetric hyperbolic system with maximal dissipative boundary conditions, well posedness of such problems is established for a large class of boundary conditions on $\partial\Sigma$. We show that our class of boundary conditions is sufficiently general to allow for a well posed formulation for different wave problems in the presence of constraints and artificial, nonreflecting boundaries, including Maxwell's equations in the Lorentz gauge and Einstein's gravitational equations in harmonic coordinates. Our results should also be useful for obtaining stable finite-difference discretizations for such problems.Comment: 22 pages, no figure

Well-posed initial-boundary value problem for the harmonic Einstein equations using energy estimates

In recent work, we used pseudo-differential theory to establish conditions that the initial-boundary value problem for second order systems of wave equations be strongly well-posed in a generalized sense. The applications included the harmonic version of the Einstein equations. Here we show that these results can also be obtained via standard energy estimates, thus establishing strong well-posedness of the harmonic Einstein problem in the classical sense.Comment: More explanatory material and title, as will appear in the published article in Classical and Quantum Gravit

The Initial-Boundary Value Problem in General Relativity

In this article we summarize what is known about the initial-boundary value problem for general relativity and discuss present problems related to it.Comment: 11 pages, 2 figures. Contribution to a special volume for Mario Castagnino's seventy fifth birthda

Fast and Slow solutions in General Relativity: The Initialization Procedure

We apply recent results in the theory of PDE, specifically in problems with two different time scales, on Einstein's equations near their Newtonian limit. The results imply a justification to Postnewtonian approximations when initialization procedures to different orders are made on the initial data. We determine up to what order initialization is needed in order to detect the contribution to the quadrupole moment due to the slow motion of a massive body as distinct from initial data contributions to fast solutions and prove that such initialization is compatible with the constraint equations. Using the results mentioned the first Postnewtonian equations and their solutions in terms of Green functions are presented in order to indicate how to proceed in calculations with this approach.Comment: 14 pages, Late

Relativistic Lagrange Formulation

It is well-known that the equations for a simple fluid can be cast into what is called their Lagrange formulation. We introduce a notion of a generalized Lagrange formulation, which is applicable to a wide variety of systems of partial differential equations. These include numerous systems of physical interest, in particular, those for various material media in general relativity. There is proved a key theorem, to the effect that, if the original (Euler) system admits an initial-value formulation, then so does its generalized Lagrange formulation.Comment: 34 pages, no figures, accepted in J. Math. Phy

Magnetospheres of black hole-neutron star binaries

We perform force-free simulations for a neutron star orbiting a black hole, aiming at clarifying the main magnetosphere properties of such binaries towards their innermost stable circular orbits. Several configurations are explored, varying the orbital separation, the individual spins and misalignment angle among the magnetic and orbital axes. We find significant electromagnetic luminosities, $L\sim 10^{42-46} \, [B_{\rm pole}/ 10^{12}{\rm G}]^2 \, {\rm erg/s}$ (depending on the specific setting), primarily powered by the orbital kinetic energy, being about one order of magnitude higher than those expected from unipolar induction. The systems typically develop current sheets that extend to long distances following a spiral arm structure. The intense curvature of the black hole produces extreme bending on a particular set of magnetic field lines as it moves along the orbit, leading to magnetic reconnections in the vicinity of the horizon. For the most symmetric scenario (aligned cases), these reconnection events can release large-scale plasmoids that carry the majority of the Poynting fluxes. On the other hand, for misaligned cases, a larger fraction of the luminosity is instead carried outwards by large-amplitude Alfv{\'e}n waves disturbances. We estimate possible precursor electromagnetic emissions based on our numerical solutions, finding radio signals as the most promising candidates to be detectable within distances of $\lesssim 200$\,Mpc by forthcoming facilities like the Square Kilometer Array

Skylight: a new code for general-relativistic ray-tracing and radiative transfer in arbitrary spacetimes

To reproduce the observed spectra and light curves originated in the neighborhood of compact objects requires accurate relativistic ray-tracing codes. In this work, we present Skylight, a new numerical code for general-relativistic ray tracing and radiative transfer in arbitrary space-time geometries and coordinate systems. The code is capable of producing images, spectra, and light curves from astrophysical models of compact objects as seen by distant observers. We incorporate two different schemes, namely Monte Carlo radiative transfer, integrating geodesics from the astrophysical region to distant observers, and camera techniques with backwards integration from the observer to the emission region. The code is validated by successfully passing several test cases, among them: thin accretion disks and neutron star hot spot emission